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We present the discovery and timing solutions of five new pulsars by stu- 
dents involved in the Pulsar Search CoUaboratory (PSC), a NSF-funded joint 
program between the National Radio Astronomy Observatory and West Virginia 
University designed to excite and engage high-school students in Science, Tech- 
nology, Engineering, and Mathematics (STEM) and related fields. We encourage 
students to pursue STEM fields by apprenticing them within a professional sci- 
entific community doing cutting edge research, specifically by teaching them to 
search for pulsars. The students are analyzing 300 hours of drift-scan survey 
data taken with the Green Bank Telescope at 350 MHz. These data cover 2876 
square degrees of the sky. Over the course of five years, more than 700 students 
have inspected diagnostic plots through a web-based graphical interface designed 
for this project. The five pulsars discovered in the data have spin periods rang- 
ing from 3.1 ms to 4.8 s. Among the new discoveries are - PSR J1926— 1314, 
a long period, nulling pulsar; PSR J1821+0155, an isolated, partially recycled 
33-ms pulsar; and PSR J1400— 14, a millisecond pulsar in a 9.5-day orbit whose 
companion is likely a white dwarf star. 

Subject headings: pulsars:general-pulsars:individual:PSR J1400— 14, PSR J1631 — 1612, 



Pulsars, rapidly rotating highly magnetized neutron stars, are most commonly st udied 



ABSTRACT 



PSR J1821+0155, PSR J1926-1314, PSR J2136-1606 



1. 



Introduction 



at radio wavelengths. The first pulsar was discovered 
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currently over 2000 of these objects have been catalogeci^ (jManchester et al" 2005 ). Over 
30,0 00 normal pulsars ar e potentially detectable in the galaxy (i.e. those beaming towards 
us) (ILorimer et al.l l2006l ) and pulsar surveys covering larger sky area with the sensitivity 
achieved with modern hardware and telescopes are the best means of detecting a larger 
sample of the pulsar population. 

Continued pulsar searches are crucial for finding exotic ob jects, which i n clude pulsars 
with ultra-fast spins (known as millisecond pulsars or MSPs; backer et al.l ( 119821 )) or in 
very relativistic orbits, binary pulsar systems, and young high velocity pulsars. Binary sys- 
tems are interesting laboratories for new physics, whether the companion is another neutron 



star, providing tests of general relativity ( ILyne et al.ll2004l : iKramer et 



20031 ). a white dwarf star, giving insight int o white dwarf cooli ng (jPriebe et al. 



al 



2004 



Burgav et al. 



19a 



or 



a main sequence star, probing stellar winds ( IMadsen et al.ll2012l ). The discovery of young, 
high velocity pulsars can s hed light on the supernova explosions in which pulsars are formed 
( lArzoumanian et al.ll2002l ). In additio n, some pulsars exhibit i ntermittent emission, such as 



Rota t ing Radio Transients (RRATs) ( IMcLaughlin et al.ll2006l ) or extreme nulhng ( ILi et al. 



20121; lEsamdin et al.l l2005l ) which give insight into the pul sar magnetosphere and how the 



magnetosphere affects the spin-down rate of the pulsar (ILyne et al.l l2010l : iKramer et al. 
20061 ). 



Pulsar searches are also necessary for finding MSPs which can be used for gravitational 
wave detection. MSPs are old pulsars that are thought to have ga ined their rapid spins during 
a process of accretion from a companion star ( lAlpar et al.lll982l ). Because of their extremely 
high rotational precision, MSPs may be able to be used to directly detect gravitational waves. 
A stochastic gravitational wave background, single, continuous gravit ational wave sources, o r 
burst sources can be detected through pulsar timing arrays (PTAs; iFoster fc Backerlll990l ). 
The success of PTAs depends crucially on the number of MSPs in the array and the sens itivity 
of PTAs increases with more MSPs distributed throughout the sky (jjenet et al.ll2006l ). 

The 100-m Robert C. Byrd Green Bank Telescope (gbtH is one of the best telescopes 
for pulsar searching because of its large collecting area, sky coverage ability, state-of-the- 
art data acquisition hardware, and location in a radio quiet zone. The GBT 350-MHz 
Drift Scan Survey, completed during the summer of 2007 while the GBT was under track 
replacement, comprised over 1491 hours of data totaling 134 TB (hereafter Papers I and 



^For an up-to-date catalog of pulsars, see |http : //www . atnf . csiro . au/research/pulsar/psrcat| 

^The Robert C. Byrd Green Bank Telescope (GBT) is operated by the National Radio Astronomy Ob- 
servatory which is a facility of the U.S. National Science Foundation operated under cooperative agreement 
by Associated Universities, Inc. 
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II, respectively; iBoyles et al.l l2012l : iLynch et al.l l2012l ). Of these data, the Pulsar Search 



Collaboratory (PSC) received 300 hours, approximately 30 TB of data. The Pulsar Search 
Collaboratory is a partnership between the National Radio Astronomy Observatory (NRAO) 
and West Virginia University (WVU) designed to engage students in Science, Technology, 
Engineering, and Mathematics (STEM) and related fields by using information technology 
to conduct current scientific research, specifically searching for new pulsars. 

One of the limiting factors in discovering radio pulsars is the large volume of data 
and the ability to analyze the diagnostic plots quickly. The data are processed through a 
software pipeline on computing clusters, searching for both periodic and transient signals. 
The resulting diagnostic plots are visually inspected, in this case, by students of the PSC. In 
its current form, plot inspection requires a basic understanding of pulsars and many person- 
hours. While the human eye is remarkably adept at distinguishing a signal from noise, the 
next generation of pulsar data analysis may require intelligent algorithms, such as neural 



next generation oi pulsar aata analysis may require intelligent aigoriti 
nets fcatough et al.lboioh . due to the large amounts of data produced. 



In this paper, we outline the data collection and processing of the entire GET 350- 
MHz Drift Scan Survey and how the data for the PSC relate to the data acquired for the 
astronomical community in ^ We discuss the PSC program in ^ In ^ we highlight the 
scientific discoveries made by the students and discuss these results in ^ 



2. Data Collection and Processing 

The data for the GBT 350-MHz Drift Scan Survey were acquired during May though 
August of 2007 while the GBT underwent a track replacement. The telescope was parked at 
an azimuth of ^ 229° for the first half of the observations and at ^ 192° for the remainder 
of the observations and covered declinations of —8° ^ S < +38° and —21° ^ S < +38°, 
depending on azimuth (Paper I). The data were collected as the sky drifted overhead and 
divided into pseudo-pointings, where a pseudo-pointing is a continuous chunk of data ^ 140 
seconds in length, equal to the transit time through the center of the beam of a source at 
350 MHz. Each pseudo-pointing overlaps the preceding one by 50%, so each piece of data is 
processed twice. This retains sensitivity to sources which may have fallen on the edge of a 
pointing (Paper II). 

The PSC received approximately 300 hours of the GBT 350-MHz Drift Scan Survey 
data. Each hour of observations results in about 48 pseudo-pointings separated by about 30 
arcminutes in right ascension, resulting in approximately 16,500 pseudo-pointings. In total, 
the PSC data covered 2,876 square degrees of the sky. Figure [1] shows the overlap between 
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the portions of the sky searched by the PSC compared to the rest of the GBT 350-MHz Drift 
Scan Survey. 

The survey was conducted at 350 MHz because of the low radio frequenc y interferenc e 
at the telescope at this frequency, the steep spectral indices of most pulsars ( jSieberlll973l ). 
and the large 36') beam size of the GBT which allows for a greater amount of sky 
coverage than highe r frequency receivers (Paper I). The Spigot autocorrelation spectrometer 
(IKaplan et al.ll2005l ) was used to acquire the data during the observations with 50 MHz of 
bandwidth at 81.92 /xs sampling time with 2048 channels and 8 bit precision. 



Each pseudo-pointing was processed using the PRESTCo software package ( iRansom 



200ll ). PRESTO performs a series of Fourier transforms to search for periodicity, including 
harmonic summing to improve sensit ivity to narrow pulses an d allowing constant accel- 
erations for pulsars in binary orbits ( ILorimer fc Kramerl l2005l ) . In addition, single-pulse 
searches for bright individual pulses are conducted. The dispersion measure (DM) is the 
integrated column density of electrons along the line of sight to the pulsar. Because the DM 
is unknown, both types of searches are made for 9700 trial values of the DM up to 1015.5 
pc cm~^. The periodicity searches are most sensitive to normal, periodic pulses, while the 
single- pulse search is designed to detect bursty, sporadic sources like rotating radio tran- 
sients ( IMcLaughlin et al.l 120061 ). PRESTO creates the periodicity search diagnostic plots 
by folding the raw data modulo the pulse period and then optimizing over local values for 
dispersion, period, and period-derivative. The students inspect the top 30 candidates output 
from PRESTO for each beam and we use PRESTO to create single-pulse search plots in 
five different DM ranges. The details of the processing pipeline can be found in Paper II. 

After the data have been processed, the students analyze the diagnostic pulsar plots 
using a quantitative grading scheme. For each pseudo-pointing, 35 diagnostic plots (30 pe- 
riodicity search plots and five single-pulse search plots) are placed in the PSC databas^. 
The database contains approximately 1000—1500 pointings visibl e at any given tim e. Only 
one student from each school can score a given pseudo-pointing ( IRosen et al.ll2010l ). Each 
pseudo-pointing must be scored by students on five individual teams for redundancy, after 
which the pseudo-pointing is no longer visible in the database and it is submitted to the 
astronomers for consideration. This system is designed to both prevent weak pulsars from 
escaping detection while also keeping the number of potential new candidates to a man- 
ageable number. Only those pointings with top rankings by students are reviewed by the 
astronomers and considered for follow-up observations. 



" ^http : //www . cv . nrao . edu/ sransom/presto/ 
^http : //psrsearch . wvu . edu 
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Right Ascension 



Fig. 1. — The portions of the sky covered by the PSC (grey), the astronomy community (hght 
grey), and where the astronomy community and PSC overlap (black). The line represents 
the Galactic plane and the circle is the location of the Galactic anticenter. 
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3. PSC Program 

A full descriptio n of the PSC p roCTa m and its educational impact on students and STEM 



careers is detailed in iRosen et al.l ( 120101 ). The PSC is currently in its fifth year and we now 
have 76 teachers and 705 students from 18 states participating. From both a scientific and 
educational perspective, the PSC has been a success. In addition to the scientific discoveries 
outlined in this paper, the PSC pr ogram has reached significant educational goals. These 



are discussed in iRosen et al.l (120101 ) and the highlights are: the PSC is reaching low-income 
students and attracts students who are first generation college-goers; the PSC succeeds in 
building confidence in students, rapport with the scientists involved in the project, and 
greater comfort with teamwork; we see additional gains in girls as they see themselves more 
as scientists after participating in the PSC program. This result is exciting and significant 
as self-efficacy is an important predictor of success in STEM fields. 

Of the approximately 16,500 pseudo-pointings in the PSC data, the students have com- 
pletely analyzed (specifically, the pseudo-pointing has been examined by five different stu- 
dents on five different teams) 6,600 pseudo-pointings, or about 40% of the data. In the 
analyzed data, the students have discovered five new pulsars, on average one new pulsar 
per 1,320 pseudo-pointings. By comparison, the astronomical community analyzed approx- 
imately 61,200 pseudo-pointings in their drift-scan data, discovering 31 new pulsars (see 
Papers I and II), resulting in a detection rate of one new pulsar per 1,900 pseudo-pointings. 
Given the random distribution of data between the astronomers and the PSC, and the ran- 
dom pseudo-pointings through which the students have already analyzed, the detection rate 
of pulsars from the PSC students is on par with the astronomers. 



4. Results 

4.1. Known Pulsars 

The presence of known pulsars in the PSC data allows us to test the survey sensitivity 
and to quantify the students' ability to detect pulsars by comparing the number of known 
pulsars expected to be detected, based on survey area and sensitivity to the known pulsars 
fiagged by the students. To determine the number of known pulsars detectable in the PSC 
data, we calculate the survey s ensitivity for a given frequ ency (u) as a function of the phase- 



averaged limiting fiux density ( ILorimer fc Kramerl 120051 . Paper II): 



C V / mm/ '^j'^ / /-I \ 
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where /5 = 1.2 is a correction factor accounting for digitization losses flLorimer fc Kramer 



20051 ) . e accounts for the drift of the source and any degradation of the signal if the source 
does not drift precisely through the beam center (Paper II), S/Nmin is the signal-to-noise 
threshold, Tgyg (K) is the combined system and sky temperature, G (K Jy""^) is the telescope 
gain, rip is the number of summed polarizations, tint (s) is the observation length, Az/ (MHz) 
is the bandwidth, and W is the pulse width. 

Figure [2] shows the students' detection of known pulsars. Each colored point represents a 
known pulsar within the sky area covered by the PSC data. The size of the point corresponds 
to the closest approach angle to the center of the 350-MHz beam with larger points being 
closer to the center. The signal-to-noise (S/N) is shown in color where we expect the students 
to detect known pulsars with a S/N > 9. The dashed lines are sensitivity curves which 
incorporate sky temperature and DM for different distances y from the beam center. 

The black squares are known pulsars that the PSC students have detected. In the 
data that the students have analyzed, there are 39 known pulsars predicted to be detectable 
above a S/N threshold of nine. Of these, the students found 32 of them as shown by the 
colored points outlined with a black square. The seven pulsars that were undetected by 
students (colored points without black squares) were likely not detectable due to strong 
radio frequency interference, scintillation, scattering, or diminished sensitivity at the edges 
of the telescope beam. These pulsars are also not detectable when we fold the raw data with 
a known ephemeris. In addition, the students have detected two known pulsars below the 
S/N = 9 threshold, represented by the two squares without colored points, and three other 
pulsars which did not have any recorded flux measurements and are therefore not plotted in 
Figure [21 We find that the students correctly found and identified all the known, detectable 
pulsars in the data they have examined. 



4.2. New Pulsars 



Once a potential pulsar was detected by the students, we conducted a confirmation ob- 
servation at 350 MHz with the GBT. If detected at 350 MHz, we then followed with a second 
observation at 820 MHz. Because of the smaller beam size a t 820 MHz, we too k a series of 
gridding observations to determine a more precise location ( IMorris et al.ll2002l ). Gridding 
entails observing a set of locations at 820 MHz which covers the area of the original 350 
MHz beam to find the Right Ascension and Declination where the pulsar has its maximum 
intensity. Upon confirmation at 820 MHz, timing observations commenced for each new 
PSC pulsar. Timing solutions provide a more precise location of the pulsar and allow us to 
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Period (s) 



Fig. 2. — The students' detection of known pulsars. The known pulsars in the area of the sky 
covered by the PSC data are shown as colored points as a function of flux density where the 
detection thres hold is S/N > 9. T he 350-MHz fluxes are calculated by assuming a spectral 
index of -1.7 ( IMaron et al.l boooh and extrapolating from 400 MHz (or 1400 MHz if no 
400 MHz values are listed) values taken from the ATNF catalog. Note that some scatter is 
expected as the observed flux of the pulsars in the PSC data may vary due to scintillation. 
The size of the point corresponds to the closest approach angle to the center of the beam 
with larger points being closer to the center. The dashed lines are sensitivity curves based 
on Equation d] for various sky temperatures and DMs at diflFerent distances y from the beam 
center. Squares represent the known pulsars that the PSC students have detected. 
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measure and/or characterize a variety of pulsar properties, such as period, period derivative, 
proper motion, and any orbital parameters, should the pulsar be in a binary system (Paper 
II). All gridding observations were cond ucted at 820 MHz u sing the Green Bank Ultimate 



Pulsar Processing Instrument (GUPPI; iDuPlain et al.l 120081 ) with 200 MHz bandwidth at 



81.92 /is time sampling resolution. Timing observations for all the pulsars except the MSP 
were conducted at 820 MHz in the same configuration as the gridding observations. For 
the MSP, we included some timing observations at 350 MHz using GUPPI with 100 MHz 
bandwidth at 81.92 /xs time sampling resolution. 

Most of the PSC pulsar timing campaigns began with a series of closely spaced ob- 
servations followed by monthly observations. Each observation lasts approximately 10-20 
minutes. From these observations, standard pulse profiles were created by summing data 
over all the frequency channels within the 200 MHz bandwidth. If the observation was at 
350 MHz, as in the case of the MSP, we used a different template. The standard profiles 
were then used to create pulse times-of-arrival (TOAs) using PSRCHIVE@. We follow the 
procedure for the timing analysis outlined in Paper I, with the modification that we create 
three TOAs per epoch for all pulsars if possible, not just binary systems. 

To calculate the fiux, we took calibration scans before and after each observation at 
each frequency. Using the fiux calibrator 3C190 measured at both 350 MHz and 820 MHz, 
we created fiux calibrated profiles using PSRCHIVE. For those observations with oflF-pulse 
fiux due to radio frequency interference, we manually estimated the on and off pulse regions 
and used these to calculate the mean fiux. 

In the data that they have examined so far, the students have discovered five new 
pulsars, one of which is an MSP. We are conducting timing observations on all of them and 
they are shown on the P-P diagram in Figure [31 The pulsars discovered by the PSC are 
indicated with stars and those discovered in the GBT 350-MHz Drift Scan Survey are shown 
with triangles. The timing parameters of the four normal pulsars, including the root mean 
square (RMS) residual ^ the difference between when we expect a pulse to arrive and when 
it actually does - are outlined in Table [421 Because the MSP, PSR J1400— 14, is in a binary 
system, we discuss its timing and orbital parameters in §4.2. 1[ Due to incomplete timing 
measurements, we refer to the MSP with only two digit precision in the declination (see 



^http: / / psrchive.sourceforge.net 
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Fig. 3. — Periods and period derivatives of all known pulsars (dots), those discovered in the 
GBT 350-MHz Drift Scan Survey (triangles), and those discovered by the PSC (stars). The 
star representing the MSP, PSR J1400— 14, has been placed at the current upper limit for 
its period derivative, since we have not yet been able to measure this quantity precisely. The 
dashed lines represent co nstant magnetic field strength and constant characteristic age (e.g. 



Lorimer fc Kramerl 120051 ). 



Table 1. Parameters of Newly Discovered Long Period Pulsars 



Parameter 1631-1612 1821+0155 1926-1314 2136-1606 



Timing Parameters 



Right Ascension (J2000) 


16:31:52.469(17) 


18:21:38.883663(30) 


19:26:53.840(16) 


21:36:00.106(68) 


Declination (J2000) 


-16:12:52.0(1.4) 


+01:55:21.98933(78) 


-13:14:04.1(1.1) 


-16:06:08.3(3.0) 


Pulsar Period (s) 


0.677683912491(22) 


0.0337813319318837(53) 


4.86428379983(10) 


1.227235400063(78) 


Period Derivative (s/s) 


1.0961(19)xl0-i^ 


2.939(24) X 10-20 


3.64393(46) X 10-1^ 


2.10(61)xl0-i^ 


Dispersion Measure (pc cm~^) 


33.768(11) 


51.75015(12) 


40.833(26) 


18.422(38) 


Reference Epoch (MJD) 


55766.0 


55584.0 


55791.0 


55764.0 


Span of Timing Data (MJD) 


55634-56128 


55584-56128 


55150-56130 


55647-56151 


Number of TOAs 


65 


84 


61 


47 


RMS Residual (/is) 


820.636 


3.432 


1580.116 


594.838 


Reduced 


8.25 


1.15 


2.95 


3.19 


Derived Parameters 


Surface Magnetic Field (10-^^ Gauss) 


0.8721 


0.001008 


13.5 


0.1626 


Spin-down Luminosity (10^^ erg/s) . 


1.39 


0.3009 


0.125 


0.004494 


Characteristic Age (Myr) 


9.8 


18200.0 


2.12 


924.0 
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4.2.1. PSR J 1400-14 

PSR J1400— 14 is the only MSP found in the P SC portion of the dri ft-scan data to 
date. At a distance of 0.5 kpc (DM of 4.6 pc cm~^; ICordes fc Lazid 120021 ) with a flux of 



4.6 mJy at 350 MHz, it is a bright, nearby pulsar. Only four other MSPs have lower DMs 
than PSR J1400— 14. At 350 MHz, we measure TOAs with accuracies of ^1 //s in one- 
minute integrations. Because of this precision, a pulse width (measured at half maximum 
from the peak of the pulse proflle) of 0.15 ms, and a location in a part of the sky with few 
other MSPs, this pulsar looks like a good candidate to be included in a pulsar timing array. 
However, at 820 MHz, we often get weak detections and sometimes do not detect the pulsar 
at all. This is likely due to the absence of bright scentless in the higher frequency band or 
perhaps because the MSP is a very steep-spectrum object. More observations are required 
to determine whether this pulsar will be useful for gravitational wave detection. 

PSR J1400— 14 has an orbital period of 9.5 days and a projected semi-major axis of 8.4 
It-s, indicating that its companion is likely a white dwarf star with a minimum mass of 0.27 
Mq (see Table [2] for the orbital parameters). The timing measurements are not complete, 
and thus we do not yet have an accurate spin-down rate or a precise location. However, by 
fitting for the position and orbital parameters, we are able to calculate an upper limit on 
the spin-down rate, as shown in Figure [31 



4.2.2. PSR J 1821+0155 



PSR J1821+0155 is unusual in that it has a 33 ms period, a spin down rat e of 2.9x10 
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s sj, a weak magnetic field of 1.0 x 10^ G (see Table [¥^21) . and no companion. iLorimer et al. 



(|2004l) interpreted such objects as being disrupted recycled pulsars (DRPs), later defined by 
Belczynski et al.l ( 1201 01 ) as isolated pulsars having 5 < 3 x 10^^ G and P > 20 ms. 



Lorimer et al.l (120041 ) described a DRP as the first born neutron star in a binary system. 



The mass accretion from the companion decreases its period while decreasing the magnetic 
field strength. The companion then undergoes a supernova explosion, disrupting the binary 
system, releasing both the neutron stars. The result is an isolated disrupted recycled pulsar 
and based on the properties of PSR J1821+0155, we classify it as a DRP. 

The "survival prob ability" of a binary s ystem to remain intact after the second supernova 
explosion is about 10% (ILorimer et al.ll2004i ). Therefore, DRPs are expected to be roughly 10 
times as common as double neutron star (DNS) binaries in the Galaxy. How ever, observations 



i ndica te that there are comparable numbers of DRPs to DNS binaries. iBelczynski et al. 



(120101 ) investigate this discrepancy using a sample of eight of the nine known DNS binaries 



Table 2. Parameters of Newly Discovered MSP 



Parameter PSR J1400-14 



Timing Parameters 

Right Ascension (J2000) 14:00:30.2422(81) 

Declination (J2000) -14:21:03.24(28) 

Pulsar Period (s) 0.003084232817869(45) 

Period Derivative (s/s) <1.9xl0"^^ 

Dispersion Measure (pc cm~^) 4.934(20) 

Reference Epoch (MJD) 56008.0 

Span of Timing Data (MJD) 55950-56121 

Number of TOAs 28 



Binary Parameters 

Orbital Period (days) 9.5476162(11) 

Projected Semi-major Axis (It-s) . . 8.419236(15) 

Epoch of Periastron (MJD) 55955.8994695(88) 

Orbital Eccentricity < 1.2 x 10"^ 

Longitude of Periastron (deg) 

Mass Function (M©) 0.007029 

Minimum Companion Mass (M©) . 0.2696 

Reduced 2.33 



Derived Parameters 

Surface Magnetic Field (10^ Gauss) < 7.8 

Spin-down Luminosity (10^^ erg/s) < 2.6 

Characteristic Age (Gyr) >0.3 
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and the 12 known DRPqj. They find that for an intermediate natal kick of a ^ 170 km 
s~^, which results in a mean velocity of 99 and 154 km s~^ for DNS and DRP populations 
respectively, there is a small 25%) observational bias toward detecting DNS binaries. The 
observational bias in favor of detecting DNS binaries is not enough to explain the discrepancy 
between the two populations, thus necessitating further study. 

We are continuing our timing observations of PSR J1821+0155. Because its profile 
is bright and narrow (see Figure Hj) , we are able to get TOAs with about 2 /xs errors for 
5-10 minute observations at 820 MHz. Continued observations may allow us to measure a 
proper motion and estimate a spatial velocity of PSR J1821+0155. Because the diflFerence in 
velocity distributions between DNS and DRP populations lead to diflFerent detection biases 
and the populations of each are small, a spatial velocity measurement of PSR J1821+0155 
can provide a significant datapoint in further exploring the discrepancy between the predicted 
and observed population sizes. 



4.2.3. PSR J1926-13U 

The pulsar PSR J1926— 1314 is notable due to its long period 4.9 s, see Table W?2\\ 
and significant nulling fraction. Nulls are times when the radio emission ceases or is greatly 
reduced. To characterize the nulling fraction in PSR J1926— 1314, we anal yzed single pulse s 



at both 350 MHz and 820 MHz, following a methodology similar to that of iRitchingsl (119761 ). 
After initially excising any radio frequency interference, we defined on-pulse and oflF-pulse 
windows, where on-pulse and off-pulse refer to regions of pulse phase based on integrated 
intensity. The on-pulse and off-pulse windows have the same number of bins, where the 
on-pulse window encompasses the entire pulse profile with as little baseline as possible. We 
then fit a second-order polynomial to the remaining bins and subtracted this from both the 
on-pulse and off-pulse windows to fiatten the baseline. To account for scintillation, for each 
200 pulse segment, we normalized the intensity by the average intensity of those 200 pulses 
(//[/]). Any single pulses with an off-pulse intensity greater than four times the RMS of the 
off-pulse windows was removed. 

Figure [5] shows a histogram of the single-pulse intensity for both the on-pulse and off- 
pulse windows at 350 MHz (left panel) and 820 MHz (righ t panel). To calculate the nulling 



fraction, we follow the methodology in iWang et al.l ( 120071 ). where we subtract the on-pulse 



^ Belczvnski et al. (|2Q1q[ ) exclude 28 isolated pulsars with 5 < 3 x 10^^ G and P < 20 ms as they are 



believed to evolved from accretion from a low mass companion, such as a white dwarf, over a long period of 
time, resulting in the evaporation of the companion. 
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J1400-14 

3.084 ms 
4.6 mJy 




J1631-1612 

0.678 s 
0.7 mJy 

I 



J1821+0155 

0.033 s 
0.6 mJy 

I 



J1926-1314 
4.864 s 
0.3 mJy 

I 



J2136-1606 

1.227 s 
1.2 mJy 

I 



Fig. 4. — The pulse profiles for the pulsars discovered by the students in the PSC drift-scan 
survey so far. All mean flux densities are from calibrated observations at 820 MHz, except 
for PSR J1400— 14, where the mean flux density is given for 350 MHz. Assuming a spectral 
index of -1.7, the equivalent flux density for PSR J1400-14 at 820 MHz is 1.1 mJy. The 
integration times for PSRs J1631-1612, J1821+0155, J1926-1314, and J2136-1606 are 
3.90, 5.89, 5.01 and 2.84 hours, respectively. The integration time for PSR J1400-14 is 15 
minutes due to the lack of observations at this frequency. The error bars on the left-hand 
side of each profile show the total expected DM smearing. 
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distribution from the off-pulse distribution, which is multiphed by a trial nulling fraction 
starting with one and decreasing incrementally until the sum of the difference of the bins 
with intensities less than zero is equal to zero. The error in the nulling fraction is the square 
root of the number of null pulses divided by the total number of pulses. The observations 
at 820 MHz show a slightly higher nulling fraction than those at 350 MHz. 



Biggsl (1 19921 ) analyzed 72 radio pulsars and found that pulsar perio d is directly propor 



tional to nulling fraction, consistent with earlier work by lRitchingsl ( 119761 ). and suggested_that 



older pulsars are harder to detect as they spend more time in their null state. I Wang et al. 



(120071 ) found that nulling fraction correlates s t rongly with large characteristic age, even more 
than pulsar period, and ICordes fc ShannonI (120081 ) showed that nulling fraction decreases 
with period derivative. 

Of the 31 pulsars discovered in the GET 350-MHz Drift Scan Survey and the PSC 
survey combined, th ree pulsars show a signif icant nulling fraction. By comparison, the Parkes 
Multibeam Survey ( IManchester et al.ll200ll ) recognized significant nulling in 23 pulsars out 
of the 750 new discoveries of normal pulsars. The high percentage of pulsars with nulling 
discovered in the GET 350-MHz Drift Scan and PSC surveys is likely due to the the sky 
area coverage being away from the Galactic plane where the pulsars are preferentially older 
and more likely to show nulling behavior. In addition, the GET 350-MHz Drift Scan Survey 
and the PSC survey are much more sensitive than previous surveys, allowing us to recognize 
nulling in more pulsars. 



5. Discussion 

One of the challenges of the PSC was the ability to manage the large number of diagnos- 
tic plots and to develop an interface for the students to quantitatively score the diagnostic 
plots. A further challenge was to collect all the scores, develop a metric for a given diagnostic 
plot, pseudo-pointing, student, or school, and to make these metrics easily accessible to the 
astronomers. The PSC database was developed for this purpose. This database, which is 
freely available, is a new scheme for analyzing the diagnostic plots that is easily scalable and 
can be applied to other pulsar searches, benefiting both astronomers and students. 

An important measure of the scientific success of the PSC is whether the students are 
able to detect pulsars as effectively as astronomers. The students have found 32 of 39 pulsars 
predicted to be detectable with S/N > 9 in the data they have examined so far and have 
re-detected an additional five known pulsars that were either below the S/N threshold or 
did not have any recorded flux measurements. Judging by detection rate, the students are 
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Fig. 5. — A histogram of the single pulse intensity of the on-pulse and off-pulse windows 
at 350 MHz (left panel) and 820 MHz (right panel) for PSR J1926-1314. The x-axis is 
the normalized intensity, which accounts for scintillation effects, and is calculated by taking 
each pulse in a 200 pulse segment and subtracting the average intensity for those 200 pulses. 
The solid lines represent the on-pulse windows and the dashed lines represent the off-pulse 
windows. The on-pulse windows show a significant number of pulses with zero or negligible 
intensity indicating a large nulling fraction. 
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on equal terms to the astronomers in finding new pulsars. We expect there to be roughly 16 
known pulsars and seven new pulsars in the remainder of the PSC data. 

Of the new pulsars that the students have discovered, only one (PSR J1400— 14) is 
an MSP. Because of the 1 /iS TOA precision at 350 MHz, narrow pulse width, bright flux, 
and location in a relatively sparsely MSP-populated region of the sky, PSR J1400— 14 is 
likely to be included in a pulsar timing array. With an orbital period of 9.5 days and a 
projected semi-major axis of 8.4 light seconds, the companion is a white dwarf star. We are 
planning optical observations to identify the companion. PSR J1821+0155 is most likely a 
DRP and therefore it will likely have received a significant natal kick from its companion. 
We will be timing this pulsar to detect proper motion in the timing residuals. This pulsar 
will likely be considered in future studies of the detection rates of DRP and DNS systems. 
PSR J1926— 1314 exhibits significant nulling and is one of three pulsars of the 31 pulsars 
discovered in the GBT 350-MHz Drift Scan Survey and the PSC surveys combined that 
show a significant nulling fraction, likely because these surveys are much more sensitive than 
previous surveys. 

The PSC is continuing to operate with the fifth summer workshop having occurred in 
July 2012. We now have 76 teachers and 705 students from 18 states participating in the 
PSC. The students are currently examining the remaining PSC data. In addition, we have 
collected an additional 50 TB of data at 350 MHz and 820 MHz using GUPPI while the GBT 
is immobilized due to maintenance or inclement weather. These data will be processed in a 
similar fashion to the original PSC data, accounting for changes in elevation and observing 
frequency. These additional data will ensure that the PSC can continue through the growth 
of additional students and teachers over the coming years. Future plans include moving 
towards a distributed project in which students from many more high schools across the 
country can participate. 

The National Radio Astronomy Observatory is a facility of the National Science Foun- 
dation operated under cooperative agreement by Associated Universities, Inc. We would like 
to acknowledge the NSF ITEST program for making the PSC possible. M.A.M. and D.R.L. 
are supported by a WVEPSCOR Research Challenge Grant and the Research Corporations. 
We are grateful for the collaboration of all of the participants in the GBT 350-MHz Drift 
Scan Survey. We would like to thank P. Demorest for useful discussions. We would like to 
acknowledge the hard work and dedication of all the PSC teachers and their high schools for 
providing support to both the teachers and students in their participation in the PSC. 
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